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Abstract 



Nonlinear cr-model is an ubiquitous model. In this paper, the 0(N) model 
where the iV-component spin is a unit vector, S 2 = 1, is considered. The 
stability of this model with respect to gradient operators (<9 M S • d u S) s , where 
the degree s is arbitrary, is discussed. Explicit two-loop calculations within 
the scheme of e-expansion, where e = (d — 2), leads to the surprising result 
that these operators are relevant. In fact, the relevancy increases with the 
degree s. We argue that this phenomenon in the 0(iV)-niodel actually reflects 
the failure of the perturbative analysis, that is, the (2 + e)-expansion. It is 
likely that it is necessary to take into account non-perturbative effects if one 
wants to describe the phase transition of the Heisenberg model within the 
context of the non-linear cr-model. Thus, uncritical use of the (2+e)-expansion 
may be misleading, especially for those cases for which there are not many 
independent checks. 
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I. INTRODUCTION 



The nonlinear a-model is an ubiquitous model that can describe systems ranging from 
quantum spins |]J to disordered electronic systems @ . In the present paper we shall discuss 
this model with O(N) symmetry. A particularly attractive method to solve this model is the 
(2 + e)-expansion ||, where the spatial dimensionality, d, defines e = d — 2. In this method, 
one takes advantage of the proximity of the non-trivial fixed point to the zero temperature 
fixed point in the limit e — * 0. One then expands around the zero temperature fixed point to 
obtain information about the non-trivial fixed point. Such arguments are clearly powerful, 
especially because there are not many explicit analytical techniques to solve this model. The 
purpose of the present paper is to examine this method more critically. 

We begin by showing that, in the most natural definition of the problem, one must 
examine the role of gradient operators of degree higher than two, although, by power- 
counting, the higher order gradient terms are irrelevant. The task of this paper is to show 
that when fluctuations are taken into account, the higher order gradient operators become 
relevant, more so as the number of gradients increases. The simplest interpretation is that 
the fixed point is infinitely unstable, and the method fails. This surprising result was 
first discovered in a related model in the context of Anderson localization of an electron 
in a random potential |3J]. However, we believe that this phenomenon in the 0(iV)-model 
actually reflects the failure of the perturbative analysis, that is of the (2 + e)-expansion. To 
appreciate this conclusion, consider the 0(3) Heisenberg model for which much is known from 
high temperature series expansion, accurate Monte-Carlo calculations, and from the (4 — d)- 
expansion of the 4 field theory. None of these methods give any hint of any pathological 
behavior. Indeed, Kehrein, Wegner and Pismak have shown that in the iV-component <f) 4 
model the one-loop contributions always make the canonically irrelevant operators even more 
irrelevant 0. Thus, it is very likely that it is necessary to take into account non-perturbative 
effects if one wants to describe this phase transition within the context of the non-linear a- 
model. At the very least, the results obtained from this expansion should be taken with 
some caution, especially in those cases in which there are not many independent checks. 

To one- loop order, the anomalous dimensions of the high gradient operators of the O(N) 
model were first calculated by Wegner ||. In an earlier paper 0, we briefly reported 
the corresponding two-loop results, and we showed that the relevance of the high-gradient 
operators persists to two-loop order. Since the method and the unexpected results may 
not be familiar, we have decided to give a full account in the present paper. We also 
analyze the renormalization group equations, and we discuss the physical consequences of 
the renormalization group flows. 

The plan of the paper is as follows. In Sec. II, we define the high gradient operators, 
and, in Sec. Ill, we set up the background field method. In Sec. IV, we calculate the 
one-loop correction to the anomalous dimension to illustrate the efficacy of the background 
field method. In Sec. V, we present the two-loop calculation. In Sec. VI, we discuss the 
flows of the renormalization group equations. Sec. VII contains our conclusions, and there 
are four Appendices. 



2 



II. HIGH GRADIENT OPERATORS 



It is useful to begin with a soft-spin model because, in many instances, the non-linear 
cx-model arises from a microscopic situation in which both the direction and the magnitude 
of the order parameter field is allowed to fluctuate. One then argues that as long as the 
O(N) symmetry is spontaneously broken, the interactions between the goldstone modes are 
precisely those given by the non- linear cr-model ||. To be specific let us consider the O(N)- 
invariant 4 field theory and explore how the high gradient operators arise. For this field 
theory, the action is given by 



(2.1) 



S[<f>] = J d d x 

Below the mean field transition, set 

<j>{x) = p(x)S(x), S 2 (x) = 1, (2.2) 

and separate out the direction, S(x), of the order parameter field 0(x) from its magnitude 
p(x). If we write 

p{x) = m + Ap(x), (2.3) 

where m is the mean field value of the magnitude of the order parameter, not containing 
any loop corrections, i.e., fluctuations. Then 

S[(j>] = J d d x + Ap(x)) 2 (d^S) 2 + X - (d^Ap) 2 + K(m + Ap(x)) 2 + u(m + Ap(x)) 1 

(2.4) 

In principle, we can integrate out the fluctuations of the magnitude of the order parameter, 
Ap(x). Because p(x) is a massive field, the effective action is local on scales larger than the 
inverse mass. It is now clear that, for scales larger than the inverse mass or the correlation 
length, the effective action will involve gradient operators to infinite order in a functional 
Taylor expansion. From symmetry, and from the locality of the action on scales larger 
than the correlation length associated with the fluctuations of the magnitude of the order 
parameter field, we find that the effective action can be cast in the form 

SeS = ^rI ddx {^ 2 + ^(^S) 4 + V^S ■ d u S) 2 + ■■■} (2.5) 

where T, f/ 4 , V4, . . . , are the coupling constants; we are measuring the temperature T in units 
of the bare spin stiffness constant. 

From power-counting, the gradient terms with powers larger than 2 are irrelevant and 
are usually dropped. If U2 S , for s > 1, denotes the coupling constant associated with the 
operator (d^S) 2s , the full scaling dimension, y 2s , of (U2 S /T) is defined by 

U2s = 2(1 — s) + e + loop corrections, (2.6) 
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where e = {d — 2); the dimension (1/T) is e, plus loop corrections. Hence, if we ignore the 
loop corrections, the contribution of the operator (<9^S) 2s vanishes as A~ 2s as the momentum 
cutoff in the theory A — > oo, for small e. In contrast, the second gradient term is relevant 
for e > and marginal for e = 0. 

The neglect of the higher order derivative terms leads to the usual partition function of 
the the nonlinear sigma model: 

Z = J DS 5(S 2 (x) - 1) exp (-± J d d x(d,S) 2 ^j . (2.7) 

This is also the naive continuum limit of the lattice Heisenberg model for which the Hamil- 
tonian is 

H/T = -- £ S, • 8 jf (2.8) 

where the temperature is measured in units of the bare spin stiffness constant. There are 
well-known arguments that the lattice Heisenberg model is in the same universality class 
as the A^-component 4 field theory. Therefore, all three models, the nonliner cx-model, 
the A^-component 4 field theory, and the lattice Heisenberg model of A^-component unit 
vector spins should belong to the same universality class, sharing the same long-distance 
and critical properties. 

However, power counting does not necessarily determine the effect of the operators. 
One must examine how loop corrections, or fluctuations, affect the picture. For the non- 
linear u model, it may appear unlikely that the corrections to the dimensions of the high 
derivative terms coming from fluctuations can overcome the canonical dimension. However, 
we shall show that, for large s, not only the one-loop correction to the canonical dimension 
is sufficiently large to render the full dimension positive, but the two loop correction is even 
larger. So, to two-loop order, the operators (<9 M S) 2s , s > 1, are relevant, contrary to the 
common description of the non-linear cr-model. 



III. THE RIEMANNIAN MANIFOLD AND THE BACKGROUND FIELD 

METHOD 



The present calculation is nearly impossible without an efficient method. In this section, 
we discuss a formalism that is efficient. The first step is to express the action as an invariant 
in the space of cosets 0(N)/0(N — 1) [[5,rU|. The length on this manifold is ds 2 = da 2 + 



{dir 1 ) 2 , but because S 2 = a 2 + tt 2 = 1, we can eliminate a using ada + 7r • dn = 0. Therefore, 
the line element is 



ds 2 = gij (it) dn l dn^ , 



(3.1) 



where the metric, gij(ir), is 



1 -7T 2 ' 



(3.2) 
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There is no unique way of choosing coordinates on this manifold. The set {ft 1 } and the set 
{n 1 } will produce two different metrics, but they are related by the transformation equation 

9l ^ ] = d^d^ 9 " 1 ^- (3 ' 3) 
The action in Eq. ( |2.7|) can now be written as 

S(tt) = ±J d d x gij (n)d^%^. (3.4) 
Similarly, the high derivative operators can be written as 

(0 w s ■ d ui s) . . . (d, s s ■ o Us s) = ( 9lJ d,yd u y) . . . ( 9pq d,yd u y) (3.5) 

Clearly, the high derivative operators and the action in Eq. ( |2.7| ) are invariant under the 
reparametrizations of the sphere, since d^n 1 transforms as a vector under reparametrization. 
The invariant measure is 



lDn]=HH,Jg(n)d7r\x), (3.6) 

i x 

where g(n) = det (&_,•). 

We briefly recall the background field method ||11|| . The strength of this method is that 
covariant expressions can be handled easily and the explicit covariance can be maintained 
at each step of the calculation. This turns to be important for an efficient organization of 
the calculations to be described. There is another important reason for using this method. 
When the fluctuations around the background field are integrated out, operators arise that 
are not invariant under the reparametrizations of the sphere. One obtains operators that 
are proportional to the classical equation of motion 

|^ = fydV + lyV^'V'': (3.7) 

On, 

where 

T ) k = ^[djgtk + d k g tj - d t g jk ] (3.8) 

are the Christoffel symbols. Since T* fc does not transform as a tensor, the equation of motion 
is not invariant, and the operators proportional to the equation of motion are not invariant 
quantities. These are the redundant operators |T2|] that do not affect the critical properties 
and can be removed by the reparametrizations of the sphere. Let us consider the shift 

n^n + 9f(n), (3.9) 

in the partition function. Here / is a smooth function, and 6 is an infinitesimal parameter. 
Then, the redundant operator, O re d, is 

t W A {'5F-£}- (3 ' 10) 
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The first term comes from the action and the second from the measure. As stated above, 
the redundant operators can be removed by a reparametrization. We adopt dimensional 
regularization and the minimal subtraction scheme for the calculation of the renormalization 
constants. In the dimensional regularization scheme, the contribution of the measure term 
can be set to zero, and the redundant operators disappear if the equation of motion is 
satisfied, which is an advantage of using the background field method. 

The anomalous dimension of a composite operator 0(x) is computed by computing the 
divergences created by the insertion in the correlation function, defined by 

= {0(x)TT(x 1 )---K(x n )). (3.11) 

The divergence of the correlation function is not only proportional to the inserted opeartor, 
but it is usually a linear combination of a number of other operators that are said to mix 
with the inserted opeartor. The renormalized operator Of is defined as 

Of = ]>>A, (3-12) 

3 

where Z^ are the renormalization constants. The renormalized correlation function Tq ^ r in 
the momentum space is given by 

r o?,fl(9»P»»0»A*) = z ^ )z ij r Oj(<l,Pi,9o), (3.13) 

where go is the bare coupling, the temperature T, and Z is the wave function renormalization 
of the field. This leads to the renormalization group equation 

r%l R (q,Pi,g,n) = 0. (3.14) 

The quantity j(g) is the anomalous dimension of the field, and 7^ is the anomalous dimen- 
sion (matrix) corresponding to the operator O, and it is defined by 

7g = -B^%(/^)^ (3-15) 

where the derivative is computed with fixed bare couplings. One then diagonalizes the 
matrix 7 at the fixed point g = g* and obtains 

r%l R (q,p i ,g;n) = 0, (3.16) 

where O a is an eigenoperator. The quantities 7 a are the eigenvalues associated with the 
eigenoperators O a , and rj = "f(g*) is the anomalous dimension of the field ir(x). The solution 
of this equation leads to the behavior 

hm r%lj q , Pi ,g*,ri = p D p- y ^ ) aR ( q /p,p l /p), (3.17) 

where D is the usual dimension of the correlation function without the operator insertion, 
and y a is given by 



d 71 
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y a = d- [O a ] - 7 a (g*), 



(3.18) 



where [O a ] is the engineering dimension of the operator O a . Of course, the largest eigenvalue 
of the matrix Jij(g*) controls the critical behavior of the operators Oj. 



In the background field method [13|, the initial step is to write the field ir(x) in terms of 
new fields: 

ir{x) = ip{x) +r](x). (3.19) 

The field ip(x) is chosen to satisfy the classical equation of motion, and 77(2;) represents quan- 
tum fluctuations. If we substitute this decomposition in the action and in the definition of 
the composite operator and expand in powers of rj, the expansion is not manifestly covari- 
ant. To obtain a covariant expansion, the field rj is written in terms of normal coordinates 
defined by the connection rf- in Eq. (|3.8| ). The normal coordinate is the tangent vector 
to the geodesic of the field manifold between if)(x) and tt(x) with a magnitude equal to the 
geodesic distance. In terms of the coordinates ip and tt, the normal coordinates, £, have the 
expansion [JT^] 

i = -n-i,+ l ~r mn {>K - ^) m (n -*/;)" + ■■■ (3.20) 

A systematic method to find the normal coordinate expansion of an arbitray operator is 
given in Ref. |13]. First, we note that for the manifold 0(N)/0(N — 1), Tfj = ir k gij, and 



the curvature tensor R™ pq is defined by the equation: 

K pq = d P K q + VT P -\p<-> <?]• (3.21) 

Therefore, 

Rmpqn gmq9np QpqQmn > (3.22) 

V 'mRpqrs = 0, (3.23) 

where V m is the usual covariant derivative on the manifold. 
For simplicity, let us introduce the notation 

GV(vr) = g mn d,7i m d u n n , (3.24) 

and we use [G ^(tt)]^ to denote the nth term of the Taylor expansion of the operator G^ iv {ti) 
in terms of the normal coordinates £. With this notation, the first four terms of the normal 
coordinate expansion are 

[GMh = 9mnWD^ m d u C + (ji~u), (3.25) 

[G^(n)} e = gmn^D^cD^ + R mi , n (^)d^ m d u ree, (3.26) 

[GWtt)]^ = ^RmiM)d^ m D^ k ee + fa <- v), (3.27) 

[G^(tt)}^ = l -R^)R lMn ^)d^ m d v r + R^D^CDuCCe; (3.28) 
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where the covariant derivative is defined by 

D,r = d,C + TT t d^C, (3.29) 

and we have used Eq. ( |3.23j ). The expansion of an invariant quantity in terms of the 
normal coordinates is manifestly reparametrization invariant: the coefficients that multiply 
the monomials £ m (x) are tensors. Note that these coefficients, the curvature tensor and the 
metric, depend on the classical field i]){x) and are not affected by the integration over the 
quantum field, £. 

The divergences that determine the renormalization constants are obtained by inte- 
grating over the quantum field £. More specifically, we calculate the one-particle irreducible 
diagrams of the expectation value 

S\d£]O ie) U £) e - s(e) &'® 

<°'"<^» = J %^m • < 330 > 

where and S ,<e - ) stand for all the terms of 0(£ n ) with n > 2 of the expansion in terms 
of the normal coordinates. The expansion of e _5<e> in powers of £ and the integration [d£\ 
generates the diagramatic expansion. 

In order to calculate the Feynmann diagrams generated from Eq. (|3.30|) , it is necessary 
to compute the propagator of the field £. However, when we use Eq. ( |3.26| ) and Eq. ( |3.29| ) 
to obtain the noninteracting part of the action, we find that 

[S} e = \f dx {g mn D,e n D»e + R lklj d mu ^d mu ^ei 1 } 

= \j dx {g mn d^ m d,C + ...}• (3.31) 

This leads to a complicated propagator that depends on the metric. Although it is possible to 
continue with the calculation, it is simpler to perform another tranformation of coordinates 
to obtain the more common propagator. The transformation is between the curved system 
of coordinates to the tangent system of coordinates: 

C = etWC, (3-32) 
C = e\e, (3.33) 



e 



"W>)g ij <* i m = S* (3.34) 



where the local matrix is known as the Vielbein. Here, we follow the convention of using 
the earlier letters of the latin alphabet (a, b, etc.) for the local indices and the latter indices 
(i, j, etc.) for the covariant indices. Furthermore, S a b is the diagonal matrix with the diagonal 
elements (1, 1, . . . , 1), so there is no distinction between covariant and contravariant indices. 
In terms of these local coordinates, the covariant derivative becomes 



= d,e + Awf^e-, (3.35) 
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where the quantity A(ip)'^ ) has dimension unity and transforms as a gauge field under the 
rotations of the tangent frames defined by e\. So, the only way in which A(%l)) ab appears in 
the calculation is through the field strength Ff£, which is 

K = d ^ - d » A t + A » A » - A t c K 

= e ai e b m i]mn d^ m d v r- (3.36) 
The substitution of the local coordinates £ a in Eq. ( |3.31| ) yields 

[S] e = \J dx {D,CD,C + elelR mtjn d^ m d,r^ b } 

= \J dx{{d,C + A a X)(d,C + A a X) + W„W^} (3-37) 
from which we obtain the usual propagator 

f p ip{x-y) 

= ^ J dpjr^, (3.38) 
In this equation m 2 is an infrared cutoff. 



IV. THE ONE LOOP CALCULATION 



We now compute the one-loop correction to demonstrate how the background field 
method works and to reproduce the results obtained by Wegner ||. A simple rescaling 
of the field £ shows that, to one-loop order, we only need to expand the action and the 
operator up to order 0(£ 2 ). Using the notation introduced in the previous section, we have 
to compute the one-particle irreducible diagrams corresponding to the expectation value 



>[S(tt)] £ 



(4.1) 



Consider the simplest case s — 2. Then 



+ [G^(n)] ( 2G aa (ip), 



(4.2) 



where the expression for [G^Jja is given in Eq. (|3.26|) and, the term [G^ir)]^ [G aa (7i)}^ 
can be found from Eq. (|3.25|) : 



(4.3) 



[G^(n)] e [G a ^7r)] e = e mc e ia d^ m d a ^D u CD a C 
+ (a «-> 0) + (/x <-> u) 
+ (a <-> f3, n <-> v). 

Note that the metric has been written in terms of the Vielbeins. Substituing in Eq. ( |4T 
and remembering that ij)(x) is a classical field, we get 
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+ ([G^)] e [G aa ^)}e)- (4-4) 

The extraction of the divergences of the one-loop diagrams is carried out with the help of 
the dimensional regularization; the details are given in Appendix A. We find 

([GM] e ) w = ^-R m n[d^ m d v r - ^d^dV], (4.5) 
([G^(7T)} e [G aP (7T)) e ) w = -^-R mpq n[^ a d^ m d^ n duVd^ q + 

(/j, <h> v) + (a <h> f3) + (/j, <r+ v, a ^ 0). (4.6) 

For the manifold 0(N)/0(N — 1), we can substitute the expression for R mnpq given earlier 
and obtain 

([G^)}e)w = l^T^W - 5 -fG aa m (4.7) 
([G, u (7r)} e [G aP (7f)} e ) w = 2G^)G a ^) - [G^G^) + (// <- v)} 

+ {a «-> {3} + {fi «-> i/} + {a «-> /3, /i «-> i/}. (4.8) 

These equations are identical to those obtained by Wegner ||. Note that the operator 
G 2 ^ u {ji) mixes with the operator G^Gp v . Furthermore, it also follows that the operators 
that mix with G mi/1 (7r) . . . G^ sUs {n) are the cyclic products of the form 

O c * c = (gijd^dpn*) (g mn dp7c m d p 7c n ) . . . {g kl d^ k d a 7i l ) (4.9) 

This is expected since within the dimensional regularization scheme, the operators that mix 
have the same symmetry and the same dimension. Surprisingly, the operators 

( fti (^ 1 )-7r i (^ 2 )^^)(^(^ 3 ) W3 7r fe (^) m4 ^) ■ ■ ■ , (4-10) 

with nil + m 2 + m 3 + • • • = 2s do not mix, even though they have the same dimension, and 
are also invariant. 

We shall be studying the renormalization of the operators O cyc . Let us imagine that 
there are N such operators denoted by the set {0^ yc }. We shall denote all other operators 
by P, and let there be M such operators. Then, the renormalized operators can be expressed 
as 

N M 

of = E^A + E^> 

3=1 3=1 
N M 

^ = E f «A + EVr (4-ii) 

3=1 3=1 
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The one-loop matrix Z^ 1 for the invariant operators with 2s gradients will have the form 



C<$ 

dm 



(4.12) 



Wegner || has shown that the largest eigenvalue comes from the matrix C^ l \ and conse- 
quently we shall not concern ourselves with the matrix D^ 1 '. 

Although we know the form of the operators that mix under renormalization, it is still a 
difficult task to diagonalize the matrix C^j since we do not have a simple way of clasifying 
the operators. Near two dimensions the problem can be solved if we introduce conformal 
coordinates HI, so that 



d + = d x + id v , d_ = d x - id y . 



In terms of these new coordinates, the divergences in Eq. Q and Eq. 
simple forms: 



jN-2) 

2vre 
(N-2) 

2ire 
1, 



7T6 



(4.13) 
J|) take very 

(4.14) 
(4.15) 
(4.16) 



while all other possibilities vanish. A vanishing result simply means that the expectation 
value is free from divergences. 

Let us now introduce the following notations: 



H = g k id + 7i k d-7T l , 

A = g kl d + n k d + 7i l i 
B = g kl d-n k d-7i l . 



(4.17) 
(4.18) 
(4.19) 



We now derive the important result that the operator H s mixes only with the set of operators 

{h s , H s - 2 (AB), H s -\AB)\ • ■ ■ , H 2 (AB)I~\ (AB)i}. 
Let us denote the normal coordinate expansion of an operator by 

Q(tt) = Q(V) + 9i (V, +Q2(M + -" ( 4 - 2 °) 
and the one-loop contribution by 

(Q(vr)) 1L = (g 2 (^,0)- (4.21) 

The last equation follows because only the quadratic term contributes to one-loop order, 
and the average implies an integration over the normal coordinates. From Eqs. ( 4.14 - 4.16 ), 
it follows that 
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(a 2 (^£))iL = iAL4(V0, 



1L 



-2/ 



H 2 ^) - A{1>)B{i>) 



(4.22) 
(4.23) 
(4.24) 



where v = (N — 2) and / = l/2ire. 

The one-loop computation, using Eqs. ( f4.22|4lT24T ), shows that 



(H s ) 



XL 



-s(s - 1)/ 



H S {^) - H s - 2 A^)B(^/j) 



(4.25) 



Similarly, 



(H s - 2 AB} lL = IH S - 2 (^)A^)B(^) [(s - l)(s - 2) - 2u] - I(s - - 2)^" 4 (^) [A(V)£(V0] 

(4.26) 

Recursively, one can show that if s mixes only with the operators listed above. It is now 
easy to show that 

((AB) l H s - 2l } 1L = - IH s - 21 (^j) [Aty)B{$)\ 1 [(s - 2l)(s - 21 - 1) - 2vl\ 

+ I(s - 2l)(s -21- l)H s - 2l ' 2 (ij) [A(^j)B(ij)} l+1 . (4.27) 

From these equations we can now derive the renormalization matrix Z^(t,e). 

Using Eq. ( |4.27| ), we can read off the matrix Zij(t, e), which is a [|] x [|] matrix: 



Z«(e,t)-1 



2vre 



a ll a n 

a 22 a 2 3 

a 33 a 34 

o o o ; ; 



where 



a i+ii+i = ( s - 2 J)( S - 2 J - 1) - 2uj, 
= ~( s ~ 2j)( s - 2 J - 1), 



(4.28) 
(4.29) 



j = 0, 1, 2 . . . , [|]. Note that there is a change in sign because to calculate Z we have to 
subtract a divergence. Because this matrix is upper triangular, the diagonal elements are 
the eigenvalues. The largest eigenvalue is the element a u . Using Eq. (|3.15| ) and ( |3.18| ), we 

get 



y£ = d-2s-^\e) 



2 + e 1 + 



8(8-1) 

N -2 



(4.30) 
(4.31) 



where we have substituted the well-known one-loop value of the fixed point t* = 2-Rtjv. 
Therefore, for sufficiently large s, regardless of how small e is, an infinite number of high 
gradient operators become relevant. 
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For the two- loop calculation we shall also need the eigenvectors of the matrix C^ 1 '. They 
are compactly contained in the matrix S that diagonalizes C^ l \ The matrix S is 



S 



(1 


Xl2 


Xl3 


X14 





1 


X 2 3 










1 













1 


V ; 









(4.32) 



where 



Ou+lQi+li+2 " " " Cii+l-li+l 



(an — ai + u + i)(an — a,i + 2i+2) • • • {&U ~ &i+u+i) 



(4.33) 



It is easily verified that SC^S 1 is diagonal and the diagonal elements are the diagonal 
elements of the matrix 



V. THE TWO LOOP CALCULATION 

The two-loop calculation is more involved. We need to subtract the subdivergences from 
the two-loop diagrams and to regulate the infrared divergences. The method for subtracting 
divergences that we follow is due to Bogoliubov-Parasiuk-Hepp-Zimmerman (BPHZ). This 
method consists of subtracting directly the subdivergences from each Feynamnn diagram by 
using the forest formula of Zimmerman ||15|| . For a complete discussion of the forest formula 



and detailed examples, we refer the reader to Collins ||16|| . With respect to infrared diver- 
gences, their presence is due to the absence of a mass term in the action. Infrared divergences 
make the computation of the loop integrals ambiguous when dimensional regularization is 
used. The poles due to the ultraviolet and the infrared divergences tend to cancel each 
other, leading in some cases to a vanishing result. The easiest way to solve this problem is 
to introduce an infrared cutoff in every propagator, i.e., we let p 2 — > p 2 + m 2 in the internal 
lines of a diagram. Of course, other choices of cutoff should not affect the final result since 
the dependence on the cutoff appears only in terms corresponding to subdivergences. These 
terms are eventually subtracted with the forest formula. 

As with the one-loop calculation, we begin by considering the normal coordinate expan- 
sion of the operator (G^Y and of the action. From a simple rescaling of the fields we learn 
that the expansion in terms of the normal coordinate needs to be carried out to order 0(£ 4 ). 
This time, a calculation of the divergences for the case s = 4 is sufficient to determine the 
divergences for arbitrary s. The different possibilities that arise in the expansion of (G^) 
are: 

0(e) : [G,M]e , [GA*)]e[G*f>{«)]& (s.i) 

0(e) :[G, '-i '-i - ' '-i 



V'l 



0(e) : [G t 



[G>i/(tt)]£3 , [G^(7r)] 5 2[G a/3 (7r)] 5 i , 
[G, u (n)] e [G aP (n)] e [G vp ] e - (5.2) 

[G^(7i)] e [G a p(TT)} e [G vp ]^ , [G Ati/ (7r)] s i[G Q! ^(7r)] $ i[G w ] $ i[G 7? p] e i. (5.3) 
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Except for the possibility of an arbitrary number of two-point insertions in the internal lines, 
Fig. [I] shows schematically the two types of diagrams that arise in the two-loop calculation. 





(b) 

w 

FIG. 1. The general form of the two- loop diagrams. The wiggly lines represent functions 
that depend on the background field tp, which is a classical field. Diagrams of type (a) lead to 
contributions proportional to double poles in e and are therefore not necessary for the computation 
of anomalous dimension. Only graphs of type (b) need to be considered. 

In Appendix B it is shown that the graphs of Fig. |l](a) are either finite or lead to double 
poles in e after the forest formula is applied to them. Hence, this type of graph plays no 
role in determining critical exponents. Moreover, the term of order 0(£ 4 ) appearing in the 
expansion of (G^Y , and in the action, can be discarded because the two-loop diagrams 
generated by the 4-point vertices are of the type shown in Fig. |l](a). In view of this result, 
the two-loop calculation requires only to carry out the expansion of an operator Q(tt) in 
terms of normal coordinates to order 0(£ 3 ), 



Q(tt) = Q(V) + 9i (V, + 92(V>, + <ki?l> , + ■ ■ 
and the only expectation values needed are: 

G(£ 3 ) '■ ([Gfiu}^[G a /3}^), ([G^,]^.[G a p]^[G K x]^) 



(5.4) 



(5.5) 
(5.6) 



In principle, it is also necessary to calculate the divergences of the operators in Eq. fl4.10| ) 
as they may mix with the operators we are considering. However, we shall show that this 
mixing may be neglected without affecting our conclusions. A detailed calculation of the 
expectation value {[G ^ lV1 ]^[G ^ 2U2 }^) is given in Appendix B. 

We emphasize again that we shall study only the renormalization of the set of operators 

{H s , H s ~\AB), H s -\AB) 2 , • • • , H 2 (AB)i-\ (AB)*} . 

The calculation is tedious, but the essential pieces are given in Appendix D. We find that 



(H s ) 



2L 



UJ 



3us + 



s(s-r 



+ 



H a - 2 {$)A{$)B($) 



S[S 



v + Q) + 2s(s -l)(s -2) 
1 



-(3-2z/) + 2s(s-l)(s-2) 



(5.7) 



where u = 24n 2 e. Similarly, 
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(H S ~ 2 AB) 2L = [9 _ v + i4( s _ 2 )] 



u 

rs-2 



+ 



if s " 2 (V')A(V') J B(^) 



(7i/-9) - (s-2)(9v+14) 



-(s - 2)(s - 3)(z/ + 2) - 2(s - 2)(s - 3)(s - 4) 



+ 



_( a _2)( a -3)(2-— ; 



+ 2(s-2)(s-3)(s-4) 



(5.8) 



The general case is even more tedious, but the necessary ingredients are given in Appendix 
D. We find that 

+ r' ll+l H s - 2l ~ 2 (il))(A(i))B(il))) l+1 + 0(4) + other operators, (5.9) 



where 



2(s - 21) (s -21- l)(s - 21 - 2) + 14/ 2 (s - 2/) 



+ 2l(s - 2l)(s - 21 - 1) + 4Z 2 (/ - 1) + -{s - 2l)(s - 21 - l){v + 6) + 6z//(s - 2/) 



-/ 2 (z/-9)-6z//-3z//(/-l) 
1 



(5.10) 



U+i 



2(s - 2l)(s -21- l)(s - 21 - 2) - 2l(s - 21) (s -21- 1), 
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+ -{s-2l){s-2l-l){s--u) 



Li-i 



UJ 



Z 2 (9-z/)+4/ 2 (s-2/)+4/ 2 (/-l) 



The matrix 7 in Eq. (|3.15|) can now be written in the form 

7 (t >e ) = -at - 26t 2 + 0(t 3 ) 

where the matrices a and b are 

/ c« \ , _ / X \ 
a - { )' { F D(2) / ' 



(5.11) 
(5.12) 

(5.13) 
(5.14) 



The matrix was obtained in the previous section. From Eqs. (|5. 10 - 5. 12 ) we find that 
the matrix is 
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12vr 2 



rn r 12 

r 2 i r 22 r 23 

r 32 r 33 r 34 

r 43 r 44 r 45 



(5.15) 



where 



—cur 



— u;r 



-uor 



(5.16) 
(5.17) 
(5.18) 



and j = 0, 1, ...[§], v = N-2. 

To calculate the anomalous dimension correct to two-loop order, it is sufficient to use 
elementary first order nondegenerate perturbation theory since the eigenvalues of are 
nondegenarate. Let U be the transformation matrix that diagonalizes the one loop matrix, 
where 



U 



The matrix S was given in the previous section. Therefore, 




U-fU 



-l 







^ o td^t- 1 



-2V 



( SC&S- 1 SXT- 1 
\ S _1 Y^T TD^T" 1 



(5.19) 



(5.20) 



The two-loop correction to the anomalous dimension y 2s is contained in the diagonal entries 
of the second matrix. This is the reason why we could ignore calculating the mixing matrices 
X and Y. We of course do not have a rigorous proof that a level crossing does not occur and 
that we do not need the block D^ 2 \ However, this is highly unlikely within perturbation 
theory. In any case, we shall show that by ignoring this block we already obtain an eigenvalue 
corresponding to a positive full dimension y 2s . A larger eigenvalue can only make things 
worse, while a smaller eigenvalue does not change our conclusions. It is easy to see that the 
we need only the (ll)-element, W-yy, of the matrix SC^S -1 , where 



Wu = [rn + x 12 r 21 ] , 



(5.21) 



with the matrix elements given above. 

Using the fixed point t — t* = ^(1 — ^), and choosing s sufficiently large, we find that, 
to two-loop order, the full dimension y 2s is given by 



V2s 



d-2s + 



es 



N-2 



l + Ol 



+ 



eV 



(N-2y 



3 V 



+ 0(e 3 ). 



(5.22) 



This is the central result of our paper. It shows that for any e, however small, we can always 
find an infinite number of high gradient operators that have positive scale dimension. 
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VI. RENORMALIZATION GROUP FLOWS 



The main result of this paper is Eq. ( |5.22 ). For any e, however small, there are an 



infinite number of high gradient operators with positive anomalous dimension y 2s . In fact, 
the dimension is larger, larger the power of the gradient operator. The two-loop contribution 
has not changed the picture obtained from the one-loop calculation 0, but has compounded 
the problem because the two-loop contribution is even larger than the one-loop contribution 
for sufficiently large s. 

The most curious phenomenon is the lack of feedback of the high gradient operators to 
the gradient operators of lower powers. It has been shown from a perturbative argument || 
that the gradient operators of power 2s contribute to the renormalization of the operators 
of powers (4s — 2). Thus, the situation is very different form the 4 theory around four 
dimensions. In that instance, the gaussian fixed point becomes unstable below 4 dimensions 
because the operator w0 4 becomes relevant. However, a new non-trivial fixed point can 
be found because of the feedback of this operator to the renormalization of the coupling 
associated with 2 term. Thus, in contrast to 4 theory, it is not possible to locate a new 
stable fixed point within the (2 + e)-expansion as described in the present paper. 

For a more complete understanding, consider the renormalization group equations. From 
Eq. ( |5.13|) , it is simple to see that 

| = _ (d _ 2)i + (W _ 2 )^ + (W _ 2) Jl_ + ... (6.!) 



^•-,- 23+ M^i) + 'v+y)] + .., (6. 2) 



U 2s dl 2tt 12tt 2 

where e l is the rescaling factor, and U 2s is the coupling constant associated with an eigen- 
operator. The equation for t is the well-known Q equation. To recover the previous result 
for y 2s , it is only necessary to substitute the fixed point value of t in the equation for U 2s . 
Written in this form, it is clear that the difficulties persist for d = 2; substitution of e = 
in Eq. ( 5.22p misleadingly leads one to believe that there are no difficulties in d = 2. The 



renormalization group flows associated with these equations for d = 2 and d = 2 + e are 
shown in Figs. ^| and || 

Let us first consider the case d = 2 + e. In principle, for t > t*, there are two possibilities: 
(1) U 2s — > 0, and (2) U 2s — > oo. The first possibility does not directly follow from these 
equations. This is the hypothetical case in which the higher order terms, not calculated 
here, bend the flows back. Because the flows are vertical as they approach t* from below, the 
magnetization must drop discontinuously at the transition |17| . Note that due to the growth 
of the couplings U 2s the spins cannot fluctuate from the preferred direction because U 2s 
multiplies the corresponding high gradient operator. Therefore, it is energetically infinitely 
costly to allow the gradient of the spins to be non-vanishing. In the second possibility, in 
which U 2s — > oo, an infinite number of U 2s grows simultaneously. This makes it exceedingly 
costly for the spins in the system to point in different directions, regardless of t. In other 
words, the spins cannot disorder. Non-perturbative effects, not contained, in the (2 + e)- 
expansion are necessary to achieve a order- disorder transition 

The case d = 2 is not much different. The growth of U 2s as t increases leads us, once 
again, to one of the previous conclusions, that is, the system cannot disorder. However, 
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t 

FIG. 2. The renormalization flows for d = 2 + e. U<is is the charge associated with the high 
gradient operator for sufficiently large s 

there appears to be a region in which the Vis initially decreases with increasing length 
scale. Therefore, for shorter length scales, the perturbative renormalization group equations 
may be approximately valid. In this region, it may be possible to match the solution 
of these low temperature renormalization group equations to strong coupling calculations. 
Effects not contained within the scheme of the perturbative renormalization group, such as 
the instanton effects described by Belavin and Polyakov ]T3 may be essential to reach a 
satisfactory description of this model. 



VII. CONCLUSIONS 



The results derived in this paper are perplexing. While (2 + e)-expansion has never 
been successful |19j in deriving the critical properties of the Heisenberg model, it has been 



useful in providing a conceptual framework. In contrast, the expansion in (4 — d) of the 
4 -fleld theory has been both quantitatively and conceptually useful. To us, the results 
derived raise serious doubts about the usefulness of the (2 + e)-expansion, because it is 
virtually certain that the phase transition in the Heisenberg model for d > 2 is described by 
two relevant operators, the temperature and the magnetic field and not by infinitely many 
relevant operators. This conviction is supported by precise finite size scaling analysis as well 
as theoretical work based on the expansion around d = 4. 

It might be argued that because the high gradient operators do not feed into the equation 
for the temperature, the disordering transition is well-described by the conventional analysis 
20fl . This argument has little force as an infinite number of parameters must be fine tuned 



to be zero. On the other hand, for small e, the pathological effect of the high gradient 
operators will be felt at very long length scales. Thus, the conventional analysis may be 
approximately valid for shorter length scales. More precise statements are difficult to make. 
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U 2s 



t 

FIG. 3. The renormalization flows for d = 2. XJis is the charge associated with the high gradient 
operator for sufficiently large s 



It is also possible that higher loop corrections may make the anomalous dimension of the 
high gradient operators negative. While this cannot be ruled out, such a situation will still 
imply rather unusual properties of the (2 + e)-expansion, if it is necessary to go to very high 
orders to eliminate the pathological behavior of this expansion. 

Based on our present understnding of the O(N) model for which there are many precise 
theoretical checks, we are forced to end on a negative note. In the theory of Anderson 
localization, (2 + e)-expansion has been used to derive a number of interesting conclusions 
concerning the distribution of the fluctuation of the moments of the conductance ET| . In fact, 



it is precisely the context in which the anomalous behavior of the high gradient operators 
|U was first discovered. While these effects may exist, it is difficult to accept them on the 
basis of the (2 + e)-expansion. The most recent work of Dupre ||22|| , in which a numerical 
simulation of a hyperbolic superplane model is carried out, indicates that there is only one 
relevant operator at the Anderson localization transition, as was conceived originally by 
Abrahams et al. 



In conclusion, we note that the (2 + e)-expansion does not seem to reproduce what 
we believe to be the correct behavior of the lattice Heisenberg model. It is possible that 
non-perturbative effects involving topological excitations are important. Such considera- 
tions have been discussed in the past [[24, HJ and have been emphasized recently in careful 
numerical simulations 125 
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APPENDIX A: 



In this appendix we provide an account of the calculation of Eqs. ( |4.5| - |4.(j| ). For conve- 
nience, we define the dimension of the fields as follows 



[ae] = K] = o, [A, 



1. 



(Al) 



Thus, any average of the field £(x) will have dimension zero. 

The one- loop divergences of the operator G^ u {ti) are obtained from the normal coordinate 
expansion through second order in the field £(x) (T^|. Making use of Eqs. ( |3.25| - |3.26| ) the 
normal coordinate expansion for the action and the operator are 



dx 



n /-a /-b 



+ A a :o{d u c + At e) + R mab nd^ m d„rei b 



(A2) 
(A3) 



Since the operator G^ u (ir) has dimension do = 2, it can only lead to divergent terms of the 
same dimension. Then, by inspection, we find that the one-loop divergence of the operator 
G^(tx) is given by 



{[G^{n)} e ) lL = D 1 + D 2 + D 3 + £> 4 
where the divergent contractions are: 

D l = [R maia2n d^ m d u r + A^A c ^}(r(x)C 2 (x)) , 



D 7 



: [RmaraM m d a r + A^ A?* ] (£ ai (x)?» {x)d^{y)d^ c {y))\ 



d 3 = -A a ;Af(c(x)d a e(x)e(y)d„e(y)) + , 

D 4 = ±A$*A?*(p(x)d^(x)e*(y)dpe'(y)d^{z)d v ?(z)) . 



(A4) 

(A5) 
(A6) 
(A7) 
(A8) 



The diagrams associated with these contractions are shown in Fig. |4j. 

Note that the divergences of these diagrams are logarithmic, so we can safely set the 
external momenta to zero. Direct and straighforward calculations give the following results: 



D 1 = I[R nm d ll 1> m d u il> n + A™A?] , 
D 2 = ~/^M a f9 a f + A%A 



D 3 
D 



-2IAfAl b , 



, ^ -I[2A ab Af + 5»„AfA° c 



(A9) 
(A10) 
(All) 
(A12) 



where I = —tt-- Thus 



2ne ' 



{[GMxl = [R mn d^ m d u r 



(A13) 
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D, 



D, 



A(¥) 





FIG. 4. One-loop diagrams contributing to the expectation value ([G^]^}. The wiggly lines 
correspond to functions of the background field. Slashed internal lines represent derivatives of the 
quantum field ^. 



It should be noted that the terms that depend on cancell out. This is a consequence 
that the only possible covariant term, F^e am e bn = Rmnijd^d^ , is antisymmetric under 
the interchange between /i and u, and then, it can not appear in the right hand side of the 
foregoing equation. 

Turning next to the derivation of Eq. (|4.6| ), we use Eq. ( |3.25| ) to obtain 



(ac <-> A) + (fJL <-> v) , 



(A14) 



where D^ a is given in Eq. ( p. 35 ). Using the fact that the dimension of the operator G^Gnx 
is d[o] = 4, we can find that the average of the first term of Eq. ( |A14| ) is given by 



(et b e na d x ^d u rD^ a D^ b )iL = E 1 + E 2 + E 3 + E, + E 5 



where 



Ex = e^d^d^A^A":^) , 

E 2 = -\e tb e na d x ^d^ n [i^A^W + A^A™*} ^{x)C\x)d il i a {y)d tl i\y)) 
E 3 = -e tb e n M%rA^A c J(^(x)d a ^x^ d (y)d K e(y)) + (^^K,a^b) , 
E 4 = -e tb e na d x ^%rAi 2b H^x)d,e(x)d^(y)d K e(y)) , 



E 5 = -e tC2 e nc A^^^ . 



(A15) 

(A16) 
(A17) 

(A18) 
(A19) 

(A20) 



The diagrams corresponding to these contraction are shown in Fig. |5|. With the exception 
of E4, all the contractions have dimension 4 and are logarithmically divergent, so we can 
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eeA(y) (3 \|/) 2 




e e (3 \\i ) 



FIG. 5. One-loop diagrams contributing to the expectation value ([G^J^i [G a p]^i). 



set the external momenta to zero. In the case of E4, the divergence is linear and, thus, 
proportional to the external momenta. In the Fourier space, the calculation of E4 yields 



q 2 (q + p) 2 



87T6 



q 2 (q + pf 



Direct calculations give the following results: 
E x = e tb e n M%rA^I , 



E- 



1 



2 = --5^Ie tb e na d x ^d u r[R p ab q d a ^ p d a r + A c a a A c a b ] , 
E 3 = 2Ie th e na AfA db d^d v r , 



£4 
E 5 



--Ie tb e na d x ^d u r(d K Af - d,A ab ) , 
~Ie tC2 e nc M%r[AfA b ^ + AfA b ? + A^Af] 



Hence, 



1 



(A21) 



(A22) 



(A23) 



he tb e na d x i;%r[(d,Af-d K Af) + (A; c Af-A% c Af^ . (A24) 

Then, we find that the terms that depend on the gauge field A ab lead to the tensor F£ b . In 
fact, it is not necessary to calculate E 5 since from gauge invariance and the result for £4, 
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the coefficient that multiplies can be known. Terms that are not gauge invariant either 
cancel or are not divergent. Making use of the definition of in Eq. ( p. 36 ), we obtain the 
one-loop divergence 

([G, u (7i)] e [G KX (n)] e ) = \ld x ^d v r\-^ K R v t nq d a Vd^ q + Rtsmnd^ m d R ^ n ] + 

(// <-> u) + (k <-> A) + (//•«-> v, k A) . (A25) 

It may be noted that the results of this appendix are valid for arbitrary Riemannian 
manifolds. 



APPENDIX B: 

In this appendix we discuss the calculation of the two-loop diagrams. First we study the 
contribution of the diagrams of the type (a) shown in Fig. [l|. Let G denote the graph in 
question, and g\ and g 2 the corresponding subgraphs. Then, if Iq denotes the Feynmann 
integral associated with G, we have that 

I G = I 9 J<n, (Bl) 
Since gi are simply one-loop graphs, dimensional regularization tell us that 

I 91 = j + k, (B2) 

where and 6j are just constants that may depend on the external momenta of the graphs 
gi. Note that, if the one-loop integral is finite, = 0. 

If T denotes the operation that selects the poles in e, the application of the forest formula 
to G yields |1J 

R(G) = I G - T(I gi )I g2 - T{I g2 )I gi , (B3) 

where R{G) denotes the overall divergence of the graph G. When I gi and I g2 are both 
divergent, the above formula yields 

R(G) = C~ + + h) - +b 2 )-^ + bl ) 

e e e e e e 

= -^ + finite (B4) 
If I g2 is finite, we obtain instead the result 

R(G) = (- + 61)62 - -62 
e e 

= finite (B5) 

From these results, we conclude that the overall divergence of the two-loop graph in Fig [l|(a) 
is proportional to double poles in e. Thus, this type of graphs do not lead to corrections to 
the anomalous dimension. 
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Now we turn to the calculation of the averages in Eqs. Q5.5H5.6Q . Since the calculation 
of the averages is quite similar, we only show the calculation of the average 

[G ^ lUl G ^ 2U2 ]i2 = {[G ^ v ^i[G ^ V2 \gi) + {[G ^ U2 \^[G ^ Vj \(2) (B6) 

Making use of Eqs. ([0§|OJ) 

[G^iuAelG^e = [e na d ul ipD^ a + ^ «-> v x ] x 

[e kb e lb D^ k D^ 1 + R kcd id^ k d^ l ^ d ] (B7) 

Substituing D^ a (see Eq. ( |3.35|) ), we obtain 

1=1 

where 

Qi{M = e^d^d^Cd^du^+^i <- vx) , (B9) 

, (bio) 

AtfAtTCT'^C + (/ii - i*)] , (Bll) 

Q 4 ty, = e^d^A^A^AtfCCT 3 + Oil - vi) , (B12) 

Q 5 (V>, = e^d^d^d^Rm^n^^d^ + (/il <- Pi) , (B13) 

Q 6 (^, EE en^^^^^^^^tana^^^e + <- ^) • (B14) 

In a similar fashion, the action 5 through third order in the field £ is given by 

[S} e + [% - fdx [St + S 2 + S 3 + S 4 + S 5 ] , (B15) 

where 

S 1 (^,0 = A a Xd a C , (B16) 

StfaQ^-JgA??? , (B17) 

^(^O^^i^V^r'f 2 , (B18) 

■W, = ^.^^ffrfef 3 (B19) 

5 B (V,0 = ^Va.^r^ffT^' • (B20) 

In principle, one has to consider all the divergent two-loop diagrams that arise from the 
contraction between one of the vertices Q n and the vertices S m . However, the number of 
diagrams can be reduced by using the dimension of the operator G Ml!/1 G M2!/2 , which is 4, and 
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by calculating only those terms proportional to the derivative of the gauge field OA. Also, 



let us recall that we are only interested in the operators defined in Eq. [4.9| . In this manner 
the only terms necessary in the calculation of [G ^ VX G ^ 2V2 ]i2 are 



w x ty) = (-i)(Qs(MSi(M) > 

W 2 (ij) = (-l)(Q 5 ^,OS^,0) , 
W 3 (i>) = (-l)(Q B (1>,0S i (il>,0) , 

WM) = (Q 5 ty, OSiiMSiiM) > 
w 5 (i>) = (Q 5 (i>, O&^OW.O) • 

From Eqs. ( |B9| - [B20| ), we see that the contractions of interest are 



'fj.ua V"itfJ \"s ^M- 



(B21) 
(B22) 
(B23) 
(B24) 
(B25) 
(B26) 



(B27) 
(B28) 
(B29) 
(B30) 
(B31) 



The diagrams associated with these contractions are shown in Fig. || To illustrate the 










A 





FIG. 6. Two- loop diagrams contributing to the expectation value [G n 2V2 ]i2 
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method of calculation, we now study the contraction A. In Fourier space, we find 

+[I El (/i 2 , yU3 )5 C3a3 5c 2 6 1 5 b2ffl2 5 Ciai + I E2 (/i 2 , Ml , /i 3 )5 (bl ^ 2) ] 
+ [I E3 (fi 2 ^ u ^)5 c ^5 b ^5 c ^5 c ^ + I Ei {^^^ {bl ^ b2) ] 
+ [a 2 <-> a 3 ] + [a 3 -> a 2 ; a 3 -> ai]} + 

{ai <-> a 2 }, (B32) 
where I Ei are just the Feynmann integrals of the diagrams of Fig. |7|. For instance the graph 






FIG. 7. The diagrams which arise from the calculation of the contraction A 



I e x yields 
Iex 



dqdk- 



<la% 3 {k - g) M2 (k-l- p) M1 (k - Z), 



[(A; — Z) 2 + m 2 ] [(A; — Z — p) 2 + m 2 } [{k — q) 2 + m 2 ] [q 2 + m 2 ] 

where we have included the infrared cutoff in the propagators. By power counting, the 
integral is linearly divergent (the superficial degree of divergence is 5 = 1). Hence, by the 
Weinberg theorem [16], the overall divergence must be also linear, i.e., the subtraction of 
divergences cancels all the non-linear divergences. Thus, we can expand the propagators of 
the previous expression in terms of the external momenta Z and p, and neglect all the terms 
except for the linear term. This yields 

k/jk^kji^k q) ^QuaQa 



iI E = (2p + 41) s / dqdk 



(k 2 + m 2 ) 3 (q 2 + m 2 



-I 



dqdk 
-(1+P)m I dqdk 



k^xik q) wQ^Qo 



(k 2 + m 2 ) 2 (q 2 + m 2 ) 

k~/(k — q) ^q^qa 



(k 2 + m 2 ) 2 (q 2 + m 2 ) 



(B33) 
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Now the calculation is reduced to computing the logarithmic divergences of the integrals of 
the right hand side of the foregoing equation. One of the ways to calculate the divergence of 
this type integral is to write the most general tensor form compatible with the symmetries 
of the integral, and then to make all the possible contractions to obtain relations between 
the coefficients. In this way we obtain 

= (2p + 4Z) a [iti<J /ll7 tf ww , + u 2 (Sij, ltl2 5 7tl3 + <^ 27 £ MlAt3 )] 
+(2p + 4Z) jU3 [w 1 5 Atl7 5 At2Q + u 2 (8 fJ , im S ja + 8^8^)) 
+ (2p + 4Z) M2 [tt3<Vi7<W + u 2 {8^^ 3 8 ia + S^S^a)] 

-\~{l)j[U4 : 5 fll ^ 2 5^ 3a -\- u^(^5 fM2 ^ 3 5^ ia -\- 6 \i x i_i 3 & fi2a)\ 

\ u ^i^^a + u${8 mll3 8 ia + S^S^a)], (B34) 

where the coefficients are 

Ul = I[_6 + e/2], «2 = -e/n, u 3 = i[6 - 5e/2], (B35) 

w 4 = 8e/fi, M 5 = 4e/fi, u B = -[-12 + 13e], (B36) 

u 7 =i[12-5e], (B37) 

and ft = 15367r 2 e 2 . 

Let us now turn to the subtraction of sub divergences of the graph I El . The three sub- 
graphs of Ie 1 are shown in Fig. Since the subgraph en is finite, the application of the 
Forest formula gives 

Rex = I El - T ( I e 12 )h 13 - T(I ei3 )I ei2 (B38) 

After the subtractions are carried out, we find that the overall divergence R El is identical 
to the integral I El except for the sign of the term . The final result for [G ^ lVl G M2i , 2 ]i2 is 
then 

[G ^ lUl G ^, 2U2 }i2 = ([G^u^^lG^^^w + ([Gmvxl^lGnsuih^w 

— \E ^m 1 nim 2 ri2 i jp X yunxn^yupq] 

L- fJ mini,m 2 n2 v\ ^2^2 ' pq,nm2 u fllfJ,2 * v\ v 2 ^ aai 

+ <-> Z^i] + [// 2 <-> U 2 ] + [/Ul <-> I/l, fl 2 <-> ^2] + 

G„ )mn «SS[^^ 5 + "J > (B39) 

where 
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e ll e 12 




ei3 

FIG. 8. The subgraphs associated with Ie 1 ■ 
E minit m 2 n 2 = ~ ^[48i? niS -R n2mirn2 + 40-R ni (n 2P ) q R pq mim2 ] , (B40) 

F pq , ni n 2 = -^[IQRpsR 3 ni n 2q ~ 4R p ( nis )tR st n2 J + K <-> n 2 ], (B41) 
G OT , wn = -^[8R ps R s mnq + 16R p{ms)t R st nq ] + [m <-> n] . (B42) 



In Eqs. ( B40 - [B42 ), we on/y show the single pole contribution, and (afe) denotes symmetriza- 



tion in the indices a and b. Note that this result is only valid for manifolds which satisfy 

^ tRmnpq 0' 

APPENDIX C: 

In this appendix we write down our results. To simplify the expressions, in this appendix 
we adopt the notation 

K^-' = d ^ mid ^ m2 ---- ( C1 ) 

The expectation values of the terms of order 0(£ 2 ) give the following results: 

<[GW]p> = -^Rna^Rj^ 3 + <OW > ( C2 ) 

l\C } i\C 1 i\ \M j\.m\n\m 2 n 2 , ju r \]> n i n 2 \]>P 9 1 _i_ 

U^/Wi/iJf 1 IV/^i^K 1 / L iKi mi«i,ni2n2 ^ ^ 2 ^2 ' PQ, n i n 2 u IJ-i ui i/ 2 aal ' 

\pti i/J + [/i 2 ffc] + [/ii <-»• 2^1, yw 2 <-> fa] , (C3) 

where 
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32 

+ -3-^mi(n lfl )^ m2 („ 2 S) *] + K ^W 2 ] , (C4) 
N p q } n in2 = —q[12 R mi ( s t) ni R m2 ni + ~^Rm 1 (n 1 s)tR m2 ( n2 

32 

+ y#mi(ms)^ m2 (n 2 S)< ] + [«1 ^ « 2 ] . (C5) 

Now we write the resuts for the expectation values of the terms of 0(£ 3 ): 

([GM = ^Rna ia2 a 3 R m Ma3 *™ , (C6) 

\(~< (~< 1 _ \rp ^r.m l n 1 m 2 n 2 , rp r \T/ rtl?l2 \T/ P 9 1 

|. (J >ii/i (J >2^2.ll2 ~~ L- tJ mini,m 2 n 2 ^ ^ 2 jy 2 1" 1 pq,n 1 n 2 u jJ,ip 2 * v x v 2 * aoJ 

+ «-> I/ x ] + [// 2 <"> f 2 ] + [A*a <-> ^1, A*2 <-» ^2] + 



where 



E mi n\,min2 ~ q [48i2 nis i2 n 2 m 1 m 2 ^Rni(n 2 p)qR Pq mi77l2 ] , (C8) 
F pq , ni n 2 = ^[IQRpsR 3 nin2 q ~ 4:R p ( nis)t R St „ 2 J + [m <-> fi 2 ] , (C9) 

Gpg imn = — ^[8-R P s-R s mn9 + 16i? p ( ms ) t i? s * n(J ] + [mi ^ rii] . (CIO) 



Also, 



([C^j^]^! [G^j^j^l [G^ 3!/3 ]^i) — pi /j,2 H3 vi v 2 "3 \Um 1 n 1 ,m 2 n 2 ,m 3 n 3 + 

+ (mi <-> m 2 ; ni <-> n 2 ) + (mi <-> m 3 ; n x <-> n 3 )] + 

^1 ^2 1^3^ aa[ u PlP2 u P3 t r \ v pq,mn 1 n 2 n 3 ) 1 

+^iM3^ 2 ^ m ( n 2 <-> n 3 ) + ^ 2M3 ^ 3 ^ m (ni «-> n 3 )]} + 
+{fx< «->!><} , (Cll) 

with 

4e 



U mi ni,m 2 n 2 ,m 3 n 3 {(-^mi(nit)n 2 ^Rmi(n2t)ni)R 



q Lv--/«i(nit)n 2 "I" °- rl mi(n 2 t)ni J-fi- n 3 m 2 m 3 

+ (n 2 <-> n 3 ; m 2 <-> m 3 )}, (C12) 

— {2-Rp( nin2 ) s -R n3gm — ^Rm{n\n 2 )sR p {n 3 q + 

-4[ J R p(ni „ 3)sJ R s „ 29m + (ni «-> n 2 )] 

-[R m (n 3 s) ni + 5R m ( ni s)n 3 )R p{n2 s) J - K <-> ^2]} • (C13) 



To find the results for the manifold 0(N)/0(N — 1), we substitute the curvature tensor, 
Eq. (|3.22|) , in the previous results. This yields, 
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192e 



a \S \b mn + \]/ mn l 



1 C-L 

p 1 T/"iini,m2n2 
q ■ r mini,m2ii2*/Jifi2i'ii'2 



+ 



12e 



ft 



where 



mini,m2Ti2 



mini 



Q pq ,ki = -l09kig P q + (11 - 3N)(g kp g ql + g kq g p i). 

Similarly, the 0(£ 3 ) expectation values are 

384e(iV-2) 
U G ^k 3 / q Jm»V ' 



f/ni /nr 1 _ V" a / miraim2ri2 1 

[ (J 'uifi ( - J 'u2^2Jl2 ~~ o A tnini,m2«2 T /ii/i2i'ii'2 1 



48e 



"1^2/ 



^2 + 



q ^pgfci * aa L"/ii i^i * ^2^2 ~ ^2 ^2 jui fi J ' 



where 
X, 



mirii,m2n2 



(20iV 66) gm\ri\ 9m2U2 iX^ < 2-N)(g rn2ni g nirni + 9m\m29ri\n\) 1 



Y pq ki = (32N - 70)g k ig pq + (23 - 10N)(g kp gi q + g p ig kq ), 
Z pqk i = (26N - 82) g kl g pq + (23 - 4N)(g kp g lq + g p ig kq ), 



Finally, 



([^Wik 1 [ G ^ 2 ]e [G^zh 1 



4e 

_ _-T ^[jmin 1 m,2n2m 3 n 3 , 

q - t miriim2ri2"i3".3 ^/ii 1/1 "2 ^2 "3 ^3 



2e 



Spqmnkl ^ 



with 



T 



minim,2n2m,3n3 



~>pqmnkl 



{(K^K^tL + + ~ ^)] + 

KJM^tL + ^2^2 J + 0* - ^3)]}, 



^32,g mini g m2n2 g m3n3 -\- 144{[(y' m3n3 ((y' mi?Tl2 (y' nira2 + fi'rrain^fi'nim^)] 

+ [m 2 <-> m 3 ; n 2 <-> n 3 ] + [mi <-> m 3 ; m <-> n 3 ]} 
-54{[#„ in2 (# 

m\mo,9m2no, ~\~ 9m2m^9m\riz)\ 

+ [m 2 <-> m 3 ; n 2 <-> n 3 ] + [mi <-> m 3 ; n x <->■ n 3 ]} 

^ ) ^{9 m Z n \9m2n^9m\n2 9m\n^9m2ri\9m-i,n2^ 1 

-104g mn g pq g k i + 64g pq (g mk g nl + g nk g m i) 

+48g m n(g P k9ql + 9qk9pl) + ^9kl{gmp9nq + 9np9mq ) 
-19{[9pn(gmkgql + 9ml9qk) + 9qn{9rak9pl + 9ml9pk)} + [m <-> n]} 
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APPENDIX D: 



In this appendix we use conformal coordinates to simplify the two-loop results obtained 
in the previous appendix. Here, we set uj = 24ir 2 e and use the following notation: 

A(tt) = d + n m d + 7i n g mn , 
B(n) = d-ir m d.ir n g mn , 
H(n) = d + 7r m d^ n g mn . 

We denote the 0(£ n ) term of the normal expansion of A(n),B(n) and H(ir) by a n (i/>,£), 
b n (i/>,£) and h n {jp,^) respectively. Straighforward manipulations yield 

(a 2 ) 2L = -—A(^) , (h 2 ) 2L = -—H(^), (Dl) 

UJ UJ 

(a iai ) 2L = A^) 2 , (ai/n) 2L = — —A(if>)H(il)), (D2) 

uj Zuj 

(aih) 2L = --[H 2 ^)(v + 1) + Aty)B(if>)(u - 1)] , (D3) 
ui 

(hih.U = -^-[H 2 (4>)(6v + 2) + A($)B($)$v - 2)] . (D4) 
We now turn to the expectation values of the terms of order 0(£ 3 ). We obtain 

(%)u = -^) , (h 3 ) = -H^), (D5) 
uj uj 

(iaa] 12 ) 2L = —A(4,) 2 , { {ah] 12 ) 2L = ^A(4,)H^), (D6) 
uj zuj 

([ab] 12 ) 2L = -\H 2 {ij){Au + 9) + A(^j)B(ij)(5is - 9)], (D7) 

UJ 

([hh] 12 ) 2L = -^[H 2 (^)(8u - 3) + A(ij)B(4>)(3 - 2v)} ; (D8) 

UJ 

where [ab]i 2 = a\b 2 + a 2 b\. Finally, 

(a.hhhL = -\H\ij) - A($)Bty)H($)] , (D9) 

UJ 

(cnhMtL = -\H 2 {4,)A{4,) - A{i>) 2 B { i>)\ , (D10) 

UJ 

{a iai h) 2L = -\H 2 ^)A^) - A(ip) 2 B(ijj)] , (Dll) 

UJ 

12 

(hhhhL = -L4(V>W)#(V0 - H\i,)} . (D12) 

UJ 

All other possibilities vanish. 
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